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*vj ' namics in the instanton liquid background, the correlation function of the O"*""*" scalar 
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Oj' quantum gluons are taken into account as well. Instead of the usual zero-width ap- 

JL , proximation for the resonance, the Brite-Wigner form for the spectral function of the 

^ ' finite-width resonance is adopted. The family of the Laplace sum rules for the scalar 

i-Q . glueball in quantum chromodynamics with and without light quarks are studcd. A 

consistency between the subtracted and unsubtracted sum rules are very well justified, 
and the values of the mass, decay width, and the coupling to the corresponding current 
for the 0^^ resonance in which the glueball fraction is dominant, arc obtained. 
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1. Introduction 

Understanding tlie nature of tlie liglitest state of glueballs, the 0"^"'" scalar glueball, is 
a long-standing puzzle in QCD[1, 2]. The mass scale of this glueball is predicted to 
be within the region of 1.30 -1.75GeV by quenched Lattice QCD[3, 4, 5, 6, 7], and by 
un- quenched lattice QCD[8, 9, 10]. Up to now, there has not been clear evidence for 
the observation of a scalar glueball, while the closest scalar candidates are /o(1370), 
/o(1500) and /o(1700)[l, 11]. More nonperturbative physics is needed in the theoretical 
and phenomenological investigation on such area. 

Laplace sum rule [12] calculations of glueball properties can be based on correlation 
functions involving interpolating field, which is not as much successful in the prediction 
of scalar glueball mass as in the ones of other hadron properties in the early days with 
the inconsistency between the subtracted[13, 14, 15, 16, 17] and unsubtracted[18, 19] 
sum rules(SSR and USSRs). Instanton vacuum should be included in the scalar 
tunnel[20, 21, 22] to offer the main non-perturbative effects. 

Instantons are the localized solutions of the classical Euclidean field equations with 
finite minimized action [23]. They can be solved by constructing the self-dual or self- 
antidual field configurations classified by different topological charges. The perturbative 
theory should be carried out around these classical solutions with the average zero- 
topological charge instead of the trivial one, which is the kernel of the semiclassical 
expansion method. 

Not all the hadrons alike[20]. Instanton contributions should not be neglected at 
least in the scalar and pseudoscalar channels. Direct instanton contributions are already 
included in sum rule approaches [24, 25, 26, 27, 28, 29] based on the instanton liquid 
model of the QCD vacuum[31, 30]. The compatibility between the resultant USSRs 
and SSR for the scarlar glueball is greatly improved, but still not be very satisfying. It 
should be noticed that, in the so-called direct instanton approximation, the interactions 
between instantons and the pure quantum gluon fields are not considered, and the 
procedure is criticized by involving with the problem of double counting[25], because in 
the correlator are included both contributions coming from condensates and instantons, 
but the latter could lead to the formation of the former. 

The interactions between the classical and quantum gluon field configurations are 
always ignored[24, 25, 26, 27, 28, 29] because the interactive effects were expected to 
be small [26]. However, there is no reasonable argument before an actual calculation. 
We have found that these interactive effects are, in fact, compatible with or even 
important more than that of the condensate and the perturbative effects at least in 
the 0"'"'' channel. Moreover, including the classical, quantum and interactive effects in 
the framework of the semiclassical expansion of the instanton background, the stability 
and the consistency for the SSR and USSRs for O"*""*" scalar glueball could be arrived 
[32]. 

Motivated by the above considerations, our main purpose in this paper is to 
investigate the O"*""*" glueball in the frame work of Laplace sum rules. To avoid the 
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problem of double counting, instead of using the scheme of the mixture of the traditional 
condensates and the so-called direct instanton contribution, we are working in the 
framework of the semiclassical expansion of QCD in the instanton liquid vacuum, 
which is a well-defined self-consistent procedure for the quantum theory justified by 
the path-integral quantization formalism. For the correlation function, we include 
the contributions from the interactions between the quantum gluons and the classical 
instanton background besides the ones coming from only instantons and from only 
quantum gluons. For the spectral function, beyond the usual zero-width approximation, 
we adopt the Breit-Wigner form for the considered resonance with correct threshold 
behavior, in order to get the information of not only the mass scale but also the full 
decay width. 

2. Correlation function 

The correlation function for the scalar glueball in Euclidean space-time with the 
virtuality q^ is defined by 

n(g2) = f d'xe''^--'{n\TOs{x)OM\^), (1) 

where 1^2) is the physical vacuum, and the scalar glueball current Og of the quantum 
numbers J^*^ = 0"'""'" is given by 

O, = a,Gl,{B)G'^'^^{B), (2) 

with Us being the strong coupling constant, it is gauge-invariant, and renormalization 
invariant at one-loop-level. On spirit of the semiclassical expansion, and in order to 
maintain the 0(4)-covariance, the gluon field strength tensor G" (i?) is considered as 
a functional of the full gluon potential 5^^ = A^^ + a^j^a with A^q and a^j_a being the 
instanton fields and the corresponding quantum fluctuations. 

The theoretical expression, n^ , for the correlation function 11 may be divided 
into the following three parts 

nQCD(g2) ^ ni^«*(g2) + ni^*(g2) + nP6"t(g2), (3) 

where Q^ = g^, and 11^^^^ (Q^), nP^^^((5^), and 11^^^ (Q^) are the contributions from the 
pure instatons, the pure perturbation QCD, and the interactions between the instantons 
and the quantum gluon fields, respectively. We note that we have not included here the 
contributions from the so-called condensates because at first in a systematic semiclassical 
expansion of QCD, the non-perturbative effects are parameterized by the classical 
instanton and anti-instanton solutions of the equation of motion of QCD, and at second 
we want to avoid the double counting problem due to the fact that some condensates can 
be reproduced from the instanton contributions, and thirdly, we have checked that the 
condensates contributions are negligible in comparison with the contributions considered 
here. 
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The perturbative contribution nP^^^((5^) up to three-loop level in the chiral limit 
of QCD is already known to be 



nPert(g2)^g4ij^ 
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where /i^ is the renormalization scale in the MS dimensional regularization scheme, and 
the coefficients with the inclusion of the threshold effects are 
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for QCD with three quark flavors up to three-loop level in the chiral limit [28, 29, 33], 
and 
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for quarkless QCD up to two-loop level [34]. Both expressions for IIP (Q^) with 
and without quark loop corrections are used in our calculation for comparison. With 
the assumption that the dominant contribution to II^'^ (Q^) comes from BPST single 
instanton and anti-instanton solutions [23, 35, 36] and the multi-instanton effects are 
negligible [26], and in view of the gauge-invariance of the correlation function, one may 
choose to work in the regular gauge of the classical single instanton potential 

'^xyjv 






xoy + p'^' 



(7) 



where r^JJj, is the 't Hooft symbol, and xq and p denote the position and size 
of the instanton, respectively. The pure instaton contribution is obtained to be 
[13, 24, 25, 37, 38, 39, 40] 



n 



inst. 



■iQ') = 2'7r'np'Q'KliVQ'p), (8) 

where K2{x) is the McDonald function, n = J^ dpn{p) and p are the overall instanton 
density and the average instanton size in the random instanton background, respectively. 
It is noticed that the contribution to 11 from the interactions between instantons and 
the quantum gluon fields is of the order of the product of a^ and the overall instanton 
density n. There is no reason to get rid off this contribution in comparison with the 
perturbative contributions of the higher order a^ considered in HP . To calculate 
such contribution, our key observation is that the instanton potential A^ obeys also the 
fixed-point gauge condition 



[X 



Xo);.A;:(x-Xo)=0, (9) 

due to the anti-symmetricity of the 't Hooft symbols. As a consequence, the instanton 
potential can be expressed in terms of the corresponding field strength tensor as follows 



A°^{x-xq)= / duuF'^^[u{x - xq)\{x - xq) 
Jo 



(10) 
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and the gauge-link with respect to the instanton fields is just the unit operator, and 
thus the trace of any product of the gauge-covariant instanton field strengths at different 
points is gauge-invariant. This allows us to conclude that the remainder quantum 
corrections to the gauge-invariant correlation function, arising from the interactions 
between the instantons and the quantum gluons, is gauge-invariant as well. Therefore, 
one may choose any specific gauge in evaluating the quantum correction to 11 . 
Working in Feynman gauge, our results for 11^ is 

nint(g2) = Coa.nvr + aln [c^ + [C^iQpf + C3] HQpf + -§-1 , (11) 



where the coefficients are: 

Co = 62.62, Ci = 1533.15, C2 = 825.81, 

C3 = -496.33, C4 = -348.89. (12) 

It is remarkable to note that the fixed point Xq, which characterizes the gauge condition 
(10), disappears in the expression of 11^ , as expected from the gauge-invariance of our 
procedure. 

3. Spectral function 

Turn to construct the spectral function for the correlation function of the scalar glueball 
current. The usual lowest one resonance plus a continuum model is used to saturate the 
phenomenological spectral function, 

Imn^HE^^) = phad^^) ^ ^^^ _ sg)ininQ^D(s), (13) 

where sq is the QCD-hadron duality threshold, p (s) the spectral function for the 
lowest scalar glueball state, and the imaginary part of the correlation function Eq. (3), 
ImnQCD(s), is 



ImnQCD(^) = - Tis' L + 2ai In ^ + f 3 In^ A, - t^A 



02 



- 167rS%pV2(pv^)F2(pv^) + aln7i{C2ph - C3). (14) 

Instead of using the zero-width approximation as usual, the Brite-Wigner form for n 
resonances is adopted for p (s) 

had^„N_V" f^rniT 



P [s 






i=l 



mf + r2/4)2 + mjT 



where /f = {fl\0 s\Of^) is the coupling of the i's resonance to the glueball current (2). 
Recall the threshold behavior for p (s) 

ff -^ AqS for s -^ 0, (16) 

with the value of Aq being fixed by the low-energy theorem of QCD[13, 18, 20, 41], 
and thus independent of what an individual resonance considered. The early QCD sum 
rule approach had often used ff — )■ Aqs (with n = 1) in the whole lowest resonance 
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region, however the obtained mass scale is too low to be expected from lattice QCD 
simulations. In fact, the threshold behavior (16) is valid in the chiral limit, it may not 
be extrapolated far away. Therefore, instead of considering the couplings fi as constants 
[26], we choose the model for / as 

\qS for s < ml, 

Xoml + A,^ for s > ml, 

with Aj being some constants, so that the spectral function p (s) has the almost 
complete Breit-Wigner form with correct threshold behavior which is important to 
maintain the convergence of the integral for the spectral function of the k = —1 sum 
rule. 



f! 



4. Finite-width Laplace sum rules 



A family of Laplace sum rules with different /c-moments can be constructed from the 
Borel transformation, B, to the correlation function (3) [12] 

£^d(so, t) = C^^^iso, t) + n(0)4,-i, (18) 

where Sq is the threshold for setting on the continuum, 11(0) is come from the subtraction 
to the corresponding dispersion relation due to the degree of divergence of the correlation 
function of the scalar glueball, and 



C^^^iso,t) = tB \i-Q')''U^^^iQ 



^^^(^o,t) 



r dss'e-'/'-p^^^{ 
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dssV^/*-nQ^^(s), (19) 
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The Laplace sum rule emphasizes the contribution from the lowest hadron state 
considered, and suppresses the higher resonance contributions and the continuum 
exponentially. 

For A; = — 1, and 1, a straightforward manipulation leads to 

c3P^{t) = [-ao + (27 - 2)ai - 0.5(67^ - I27 - 71^)02] t^ + 2^^^ 
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^QCD^^^ = [_6ao + (127 - 22)ai - (187^ - 667 - Svr^ + 86)02] ^^ 



+ 2^Tx'^n—e 



(9 - Ax)Kq{x) + ( 7 - 4a; + - ) Ki{x) 



X 

2+3 



+ nai[2C2p'r - C^t\ (23) 

where x = p^t/2. 

5. Numerical calculation 

Now, we specify the input parameters in numerical calculation. We take the color and 
flavor numbers to be iVc = 3 and Nf = 3, respectively. The expressions for two-loops 
quarkless {Nf = 0) running coupling constant as{Q'^) at renormalization scale /i [42, 43] 
and for the three-loop running coupling constant with three flavors {Nf = 3) are used, 
where the central value of the MS QCD scale A is taken to be 120 MeV. We recall 
here that a research on the renormalization group improvement for Laplace sum rules 
amount to choose the renormalization scale to be /i^ = t [44]. The subtraction constant 
n(0) is determined by low-energy theorem [13] 

397r 
n(0) = ^^{asG^) ~ 0.6 GeV^ (24) 

which leads to Aq = S.OGeV. The values of the average instanton size and the overall 
instanton density are adopted from the instanton liquid model 

n = lfm-^ = 0.0016 GeV^ p = -fm = 1.689 GeV~\ (25) 

Finally, the mass of the neutral pion is taking from the experimental data, i.e. m^^ = 135 
MeV. 

To determining the values of the resonance parameters appearing in Eq.(15), we 
match both sides of sum rules (18) optimally in the flducial domain. The conditions 
for determining the value of Sq are: flrst, it should be grater than m"^; second, it should 
guarantee that there exists a sum rule window for our Laplace sum rules. We note 
that the upper limit tmax of the sum rule window is determined by requiring that the 
contribution from the continuum should be less than that of the resonance 

^^''Hsctmax) < /:?^^(^o,tmax), (26) 

while the lower limit t^^^ of the sum rule window is obtained by requiring the 

OCD 

contribution of pure instantons to be greater than 50% of C^ (so,^), because such 

classical contributions should be dominant in the low-energy region. Moreover, to 
require that the multi-instanton corrections remain negligible, we simply adopt a rough 
estimate 

^m\n ^ (2p)' ~ (2/0.6GeV)l (27) 

In order to measure the compatibility between both sides of the sum rules (18) realized 
in our numerical simulation, we introduce a variation, 6, deflned by 

1 ^ [L(t.) - R{t,)r ,.„. 
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where the interval [tj^^jj^, tmax] is divided into 100 equal small intervals, A^ = 101, and 
L{ti) and -R(tj) are l.h.s and r.h.s of Eq.(18) evaluated at tj. 

In the world of quarkless QCD there is only one well defined scalar bound state of 
gluons below IGeV suggested by lattice QCD, and thus we choose n = 1 in Eq.(15). 
Including quarks enhances the difficulty of the task since many states possessing the 
same quantum numbers may present in the correlator. The assumption of a single well- 
isolated lowest resonance is questioned from the admixture with quarkonium states, 
and from the experimental data that three O"*""*" scalar states around the mass scale of 
1500MeV, namely /o(1370), /o(1500) and /o(1710). Therefore, we choose n = 3 for 
p '^{s). With the requirements mentioned above, the optimal parameters governing 
the sum rules are listed in Tab.l. The corresponding curves for the l.h.s. and r.h.s. of 

Table 1. The fitting values of the mass m(GeV) and width r(GeV) of the 
possible 0^"'" resonances, and of the parameters A(GeV), /(GeV) and so(GeV^) 
characterizing the couplings to the corresponding resonances and the continuum 
threshold, as well as of the sum rule window of t(GeV^) and the compatibility 
measure 5 for finite- width Laplace sum rules (18) in the quarkless QCD (shown 
in the first three lines), as well as in QCD with three massless quarks (shown 
in the remainder lines, where the first three lines for n = 1, and the others for 
n = 3). 
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(18) of fc = —1, and -|-1 in quarkless QCD for n = 1, and in QCD with three massless 
quarks for n = 3 are displayed in Figs.l and 2, respectively. These figurations show the 
consistent match between the both sides of Eq. (18) for k = —1,0 and 1 respectively, 
with the corresponding fitting parameters. The solid lines are the r.h.s. (QCD) of Eq. 
(18), and the dashed lines are the l.h.s. (HAD) of Eq. (18), while the dotted line for the 
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r.h.s. (QCD) excluding the contribution of interactions between the instantons and the 
quantum gluons. The matching between both sides of the sum rules are very well over 
the whole fiducial region with a very little departure. In the case of QCD with three 
massless quarks, the curves for n = 1 are similar to those for n = 3 with little worse 
compatibility, and not displayed here. 

6. Conclusion and discussion 

The properties of O"*"^ glueball are examined in a family of the finite-width Laplace 
sum rules. The correlation function is calculated in a semiclassical expansion, a well- 
defined process justified in the path-integral quantization formalism, of QCD in the 
instanton background, namely the instanton liquid model of the QCD vacuum. Besides 
the contributions from pure gluons and instantons separately, the one arising from the 
interactions between the classical instanton fields and the quantum gluon ones are taken 
into account as well. Instead of using the usual zero- width approximation for the spectral 
function of the considered resonances, the Breit-Wigner form for the resonances with a 
correct threshold behavior is adopted. With the QCD standard input parameters, three 
Laplace sum rules with k = —1, 0, 1 are carefully studied. 

By taking the average for the values, listed in tab. 1, of the corresponding k = —1, 
and 1 sum rules, for the case oi n = 1, the values of the mass and width, and the 
other optical fit parameters are 

m = 1.51 ± 0.15GeV, F = 0.12 ± 0.04GeV, 
/ = 1.52 ± 0.12 GeV, Sq = 4.5 ± 0.5 GeV^ 

in quarkless QCD, and 

m = 1.50 ± 0.15GeV, F = 0.12 ± 0.04GeV, 
/ = 1.61 ± 0.11 GeV, So = 4.5 ± 0.5 GeV^ 

for QCD with three massless flavors, where the errors are estimated from the 
uncertainties of the spread between the individual sum rules, and by varying the 
phenomenological parameters, A and {agG"^), appropriately away from their central 
values A = 120 - 200 MeV and (agG^) = 0.6 - O.SGeV^ It is remarkable to notice that 
these two results are close in value, and it indicates that the considered quark effects 
may be not so large at the energy scale of the resonance mass above IGeV. 

The above conclusion is further justified in the investigation of the case oi n = 3. 
We can see that the current is coupled mainly to the resonance state near /o(1500) 
predicted in a single resonance approach. To be quantitative, let us consider the case 
with the most excellent compatibility, namely the k = results shown in Tab.l and 
Fig.2. The corresponding couplings /^ to the three resonances Ri, R2 and R3 with 
masses 1.35GeV, 1.47GeV and 1.70GeV are 

0.092GeV^ 4.073GeV^ 0.094GeV2 (29) 
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Figure 1. The curves for the l.h.s. and r.h.s. of the Eq. (18) for k = —1, and 
+1 cases of quarkless QCD with only the lowest resonance considered. The 
sohd Hne denotes the r.h.s. (QCD), dashed hne for l.h.s. (HAD), and dotted 
line for the r.h.s. (QCD) the without interaction contribution. 
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Figure 2. The curves for the l.h.s. and r.h.s. of the Eq. (18) for k = —1, 
and +1 cases of QCD with three massless quarks and the lowest three 
resonances considered. The solid line denotes the r.h.s. (QCD), dashed line 
for l.h.s. (HAD), and dotted line for the r.h.s. (QCD) the without interaction 
contribution. 
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respectively, Note that 

(O|0,|i?i) \ / (O|0,|l) \ 

(0|O,|i?2) =M (0|O,|8) (30) 

(o|o,|i?3) / V <oi^^l^) / 

where M stands for the mixing matrix (s. the second one of Eqs.(45) in Ref.[ll])- The 
values of the couplings of Og to the qq states |1), |8) and the pure glueball state \G) are 

(O|0,|l) = -0.48GeV^ 

(0|O,|8) = 0.19GeV^ (31) 

(0|O,|G) = 0.86GeV3 

after normalization, respectively. Although the estimation is relatively rough, it is still 
remarkable to notice that, first, the coupling to \G) is dominant; and second, the signs 
of the couplings to |1) and |8) is consistent with the scalar glueball- meson coupling 
theorems[13, 20, 51]. 

As summary, we may conclude that the values of the mass and decay width of 
the O"*""*" resonance, in which the fraction of the scalar glueball state is dominant, are 
m = 1.5 ±0.15 GeV and F = 0.12 ±0.03 GeV, respectively, and the value of its coupling 
to the corresponding current is Xs>mi = 1.5 ± 0.12 GeV. They are not only compatible 
with lattice QCD simulation [3, 4, 5, 6, 7] and other estimation[24, 25, 26, 46, 47], but 
also in good accordance with the experimental data of /o(1500) [48, 49, 50]. 

It is also remarkable that the three Laplace sum rules lead to almost the same 
results, a consistency between the subtracted and unsubtracted sum rules are very well 
justified. We note that we have not working within the mixed scheme, namely with 
including condensates, and in the same time, adopting the so-called direct instanton 
approximation, but simply with a self-consistent framework, a quantum theory in a 
classical background, without the problem of double counting. In this aspect, our results 
further justified the instanton liquid model for QCD among other many justifications. 

In our semiclassical expansion, the leading contribution to the sum rules comes 
from instantons themselves, especially in the region below the threshold sq. It is the 
amount of this contribution determining the low-bound of the sum rule window. This 
means that the non-linear configurations of gluons have a dominant role with respect 
to the quantum fluctuations in the low-energy region. 

The contribution of the interactions between the classical instanton fields and 
quantum gluon ones, considered in this letter but neglected in earlier sum rule 
calculations [24, 25, 28, 29, 46], is in fact not negligible. To the contrast, its amount 
is approximately double or even triple of that from the pure quantum fluctuations in 
the whole fiducial domain, expected from a view point of the semiclassical expansion. 
Moreover, it is obviously seen from Figs. 1-2 that, without taking the contribution from 

the interactions between instantons and quantum gluons into account, the departures 

1, ^ /'had/ J.N J ^QCD without interaction/ ,-, -, , j n .i 

between L^°'^(so,tj and L^r {so,t) become large, and all the 

three Laplace sum rules become less stable, and thus less reliable. 
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Finally, it should be noticed that the imaginary part of instanton contribution is an 
oscillating, amplificatory and nonpositive defined function, and so is the imaginary part 
of the correlation function. This property which is a fatal problem for the QCD sum rule 
calculation with the instanton background, may make the contribution of continuum 
too large to be under control. Hilmar Forkel introduced a Gaussian distribution for 
instanton to get rid off this trouble, and obtained a smaller 0+"'' mass scale: 1.25 ± 0.2 
GeV [24] compared to the earlier result 1.53 ±0.2 GeV [25]. We didn't use this Gaussian 
distribution, but simply choose a smaller fitting parameter to avoid this problem. 
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